clear;clc;

breaks      =   linspace(-1,1,20)';
evennum     =   0;
k           =   3;
fspace      =   fundef({'spli',breaks,evennum,k});
node        =   funnode(fspace);
value       =   node.^3;
order       =   0;
B           =   splibas(breaks,evennum,k,node,[0;1;2]);
Coef        =   B{1}\value;

X           =   linspace(-2,2,30)';
B_test      =   ExtendedSplineBasis(breaks,evennum,k,X,2);
plot(X,B_test*Coef);